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1 Introduction 



Starting from the seminal paper by E.Lieb, T.Schultz and D.Mattis the study 
of the one-dimensional spin-| XY models has attracted much interest. A lot 
of exact results concerning thermodynamics, spin correlations and their dy- 
namics were found over the last 35 years. The analytical results obtained 
for random versions of such models are not so impressive. Some results deal- 
ing with magnetic properties for special cases of XY models with random in- 
tersite interactions were found by E.R.Smith/ E.Barouch and B.M.McCoy^ 
and R.O.Zaitsev^ using an approach developed by F.J. Dyson. A somewhat 
different approach was suggested by H.Nishimori/^ who presented exact cal- 
culations of thermodynamic quantities for the isotropic XY model in a ran- 
dom lorentzian transverse field. Nishimori's exact solution is based on the 
Jordan- Wigner transformation of the spin Hamiltonian to a tight-binding 
model of non-interacting spinless fermions with diagonal lorentzian disor- 
der. For the latter fermionic model the random-averaged one-particle Green 
functions (and hence the density of states that yields thermodynamics) were 
found exactly first by P.Lloyd^^ with the help of contour integrals. Later on 
Nishimori's work was generalized for models with alternating bonds^^ and 
additional intersite Dzyaloshinskii-Moriya interaction.^^ The main results ob- 
tained in the papers on models with a lorentzian transverse field concern the 
changes in the temperature dependences of entropy, specific heat and static 
transverse linear susceptibility as well as the ground-state transverse mag- 
netization as a function of averaged transverse field once the randomness is 
introduced. 

The idea of the present paper is to study the thermodynamics of a ran- 
dom spin-i XY chain exploiting an extended version of the Lloyd model 
with off-diagonal disorder^^. In the late seventies in a series of papers this 
extended Lloyd model was applied to disordered electron systems by one of 
the authors. Similarly to Ref. 11 we present exact calculations for vari- 
ous thermodynamic quantities. However, due to the additional off-diagonal 
disorder our results concerning the infiuence of disorder on thermodynamic 
functions differ to some extent from those obtained by H.Nishimori^^. 

The paper is organized as follows. In Section 2 we describe the Jordan- 
Wigner transformation from spins to non-interacting spinless fermions and 
the evaluation of the averaged one-fermion Green functions. In Section 3 
the averaged Green functions are used to calculate the thermodynamic prop- 
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erties, namely, entropy, specific heat, transverse magnetization and static 
transverse linear susceptibility. Some conclusions are given in Section 4. A 
short report of these results was presented in Ref. 21. 



2 Jordan- Wigner transformation 

and averaged one-fermion Green functions 

We consider a linear XY chain of N spins | in a transverse field with peri- 
odical boundary conditions. The Hamiltonian reads 

N N 

n=l n=l 

N , 1\ ^ J 

= X! ( ^n^n ~ ) + X! 17 i^n^n+l + *n ^n+l) ) ^n+N — ^n^ (l) 
n=l ^ n=l ^ 

where Qn is the transverse field at site n and J„ is the exchange interaction 
between the sites n and n + 1. The J„ are taken to be independent random 
variables with a lorentzian probability distribution 

pW = I- ^TT-TT,- (2) 

^ [Jn - Jo) + 

Here Jq is the mean value and F is the width of the distribution (strength of 
disorder). In order to treat the model (1), (2) in an exact manner we assume 
the following relation between the transverse field at each site Qn and the 
surrounding intersite interactions (cf. Ref. 15-20) 

Q„ — Qq = o ~ "I — ~T, — ~1 ' ^ is real, | a |> 1, (3) 



where Qq is the averaged transverse field at site. 

Next we transform the spin model to a fermionic model by Jordan- Wigner 
transformation: Ci = s{ , cf = sf, Cj = Pj^isj , cf = Pj^isf, j = 2,...,N, 
Pj = n;!(=i(— 2s^). The resulting Hamiltonian reads 

H^H- + BP+ = H+P+ + H-P-, 
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^ N N N J 

^ n=l n=l n=l ^ 

B = -Jn (c%ci - cncI^ , 

lip N 

p^^^^ p^n(-2<) (4) 

^ n=l 

with anticyclic boundary conditions for and cyclic boundary conditions 
for H~ . For the calculation of thermodynamic properties of model (1) one can 
omit the boundary term B,^^ i.e. it is sufficient to study the thermodynamics 
of spinless fermions described by the c-cyclic Hamiltonian H = H~ . This 
c-cyclic fermionic Hamiltonian corresponds to the one-dimensional version 
of Anderson's model with off-diagonal disorder considered by W.John and 
J.Schreiber.^^ 

Following Ref. 15 one introduces the retarded and advanced temperature 
double-time Green functions G^^{t) = ^i9{±t) < {c„(t),c+} >, G^^(t) = 
^ J^^ doje-''^*Gl^{uj ±i£),£^ +0. For Gl^{uj ± is) one finds the following 
set of equations 

{Uj±ie- Qn) Gnmi^ ± " 

Suppose that J„, ... (and hence ...,f2„, ...) are complex variables. The 
singularities of G^^(a;±i£) are given by the zeros of det{uj±ie — H). det(a;± 
ie—H) is different from zero if the eigenvalues A of Im(t(;±z£— i?) are either all 
positive or all negative. /,From the Gershgorin criterion^^'^^ for the complex 
matrix lm(a; ±ie — H) one gets that for any eigenvalue A at least one of the 
conditions 

1 1 

I Im(a; ±ie — Q„) — A |< - | ImJ„_i | -|-- | ImJ„ | (6) 

has to be fulfilled. Using (3) the inequalities (6) can be written as 

I lm{u ±ie) (IniJ„_i + ImJ„) - A |< ^ | ImJ„_i | +^ \ ImJ„ |, | a |> 1. (7) 

Let us consider the retarded Green function {lm{uj+ie) > 0). Then according 
to (7) for a > 1 all A must be positive if all ImJ„ < 0, whereas for a < —1 
all A must be positive if all ImJ„ > 0. Similarly, for the advanced Green 



J, 



n-l 



(5) 
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function (Im(a; — ie) < 0) according to (7) for a > 1 all A are negative if all 
Im Jn > 0, and for a < —1 all A are negative if all ImJ„ < 0. 

Consequently, for a > 1 (a < 1) the retarded Green function G~^{uj + i£) 
cannot have a pole in the lower (upper) half-planes of complex variables J„, 
whereas the advanced Green function G^j^{uj — is) cannot have a pole in the 
upper (lower) half-planes of complex variables J„ for a > 1 {a < —1) . Using 
these properties one can perform the averaging of equations (5), defined by 

^ roo 1 r 



(•••)- n/ dj^-T-j — 7^^(-) 

n/ ^^n-7T-T^7FV7T-T— T^(-)> (») 

n=l 



oo TT ( Jn - Jo + iV) {Jn - Jq - iT) 



by means of contour integrals. For the averaging of a function F(..., O^, J„, ...) 
that has no poles in lower half-planes J„, one can close the contours of inte- 
gration in (8) in these half-planes. One obtains 



F(..., Q„, J„, ...) = F(..., Qo - iaV, Jo - iV, ...). (9) 

Similarly, for the function without poles in upper half-planes J„ one gets by 
contour integration 



F(..., Jn, ...) = F(..., 1^0 + Jo + tV, ...). (10) 
Then the averaged equations for Green functions (5) due to (9), (10) read 



[a; ± z£ - (Qo T ^ I a I V)] Q^(u; ± le) - 

JoTisgn(a)r 



Gt-l,m{^ ± + G^+l,m{^ ± = Snm- (H) 

Equations (11) possess translational symmetry and therefore they can be 
solved in a standard way. The resulting averaged Green functions read 

1 



=F i I o I r -|- ( Jo =F isgn(a)r) cos k] 

-(12) 



with X = u) — VLq ± i \ a \ V , y = Jq ^ isgn(a)r. 
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3 Entropy, specific heat, transverse magneti- 
zation and static transverse linear suscep- 
tibility 

The obtained averaged Green functions (12) allow to study thermodynam- 
ics of the spin model (l)-(3). For this we diagonalize the bilinear in Fermi 
operators form H~ (4) by the canonical transformation -q^ — J2n=i QknCm 

■^kdkn — ^i=l Qki-^iny ^ij ~ ^i^ij "I" 2'^«*^i>i+l + 2Ji—l^j,i-li ^i=l 9ki9pi ~ ^kp-i 

Yjp=i9pi9p3 = Sij with the result Y.p=i^piVpVp ~ |)- The thermodynam- 
ics for certain realization of random intersite interactions is then deter- 
mined by the spectrum of elementary excitations Ap or its density p{E) = 
J2p=i S{E — Ap). For example, the Helmholtz free energy per site is given 
by / = ^{-^lnnpli[exp(-^) +exp(^)]} = -1 /d£;p(£;) ln(2chf ). 
The result of averaging over realizations of random variables is given by 
/ = J dEp{E) ln(2ch^). The required averaged density of states is then 
calculated by the averaged one-particle Green functions (12) 



^— Im 



p{E) = --lmG-^{E) = -ImGUE) 



TT 



TT 



1 



TT ^/fE -no±i \ a\Tf - {JoT isgn(a)Tf 



1 

TT 



VA^ + B^-A 
\ 2(A2 + S2) ' 



A={E- flof + (1- I « r)r' - Jo, B = 2r[| a\{E-flo)+ sgn(a) Jo].(13) 



Really,p(E) = ^E 



N 



-^lmT-^{E + te) 



AT 2^p=l 



^imr+(i? 



ie) , where 



TJ(t) = TtO{±t) < {Vp{t),vt} >= Ei=l^J=l9p^9qJGfJ{t), and therefore 



p{E) = jrEf=i -ilm GjAE + te 



E-ie) 



the aver- 



N ^3=1 

aging yields the first two equalities in the left-hand side of (13). 

Knowing the averaged Helmholtz free energy we can calculate the entropy 
and specific heat by the formulae 



0- 



dp 



J dEp{E) 



In 2ch 



(5E^ 



2 2 



(14) 
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-^-^l,E,iE)[-^) . (15) 



Due to the magic property of (13) -^p{E) — —-^p{E) one can express 
transverse magnetization and static transverse linear susceptibility through 
the density of states 



1 ^ 

rn;, = < 



Let us discuss the obtained results. 

Note at first, that in the absence of randomness (13) reduces to the well- 
known result 

if \E-no\<\Jo 



p{E) ^^-Im , ={ " V^i-i^-^^r ' "(18) 

^ J{E -no±isf - ( otherwise 



as anticipated. 

The considered model (l)-(3) essentially differs from that one with diag- 
onal disorder treated by H.Nishimori^^ ( J„ = J, f2„ are independent random 
variables with lorentzian distributions). However, Nishimori's model can be 
obtained as a certain limit of the present model, namely F — > 0, | a | F = 
const = Fjs^. The density of states (13) in contrast to the case of diagonal 
disorder is not symmetric with respect to the change E — Qq ^ —{E — Qo). 
However, it remains the same after the replacement E — Qq ^ —[E — Qq), 
a — > —a, OT E — D,o ^ —{E — fio), Jq — Jo, since the simultaneous change 
of signs of Jo and a in (13) does not affect p{E). Without loss of generality we 
choose fioi Jo > throughout the rest of the paper. It is also convenient, al- 
though by no means essential, to put hereafter Jq = 1. The above-mentioned 
symmetry of the density of states can be seen in Figs. 1-3, where the av- 
eraged density of states (13) (Fig. 1), the averaged density of states (13) 
in comparison with histograms p{E) calculated for a certain realization of 
random intersite interactions using exact finite-chain calculations^^ (Fig- 2), 
and the histograms p{E) obtained by the latter approach for | a |< 1 (Fig. 3) 
are displayed. The density of states for the non-random case (18) is depicted 
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in Fig. 1 by dashed lines. For large | a \ the edges of the zone are completely 
smeared out with increasing strength of disorder F; for | a |~ 1 the increase 
of disorder results in a smearing out of mainly one edge of the zone. A further 
decrease of | a | up to leads to a recovering of the symmetry with respect 
to E — Qq -~{E — Qq) and to transforming of p{E) into 5{E — Qq). Some 
consequences induced by this dependence of p{E) on F and a will be seen in 
the behaviour of thermodynamic quantities. 

The results of numerical calculations of thermodynamic quantities are 
presented in Figs. 4-13, namely, the temperature dependences of the entropy 
(14) (Figs. 4,5), the specific heat (15) (Figs. 7,8), the transverse magnetiza- 
tion (16) (Fig. 11) and the static transverse linear susceptibility (17) (Fig. 
13), the dependence on averaged transverse field at low temperatures of the 
entropy (Fig. 6), the specific heat (Fig. 9), the transverse magnetization 
(Fig. 10) and the static transverse linear susceptibility (Fig. 12); the curves 
that correspond to the non-random case are depicted in Figs. 4-13 by dashed 
lines. 

The influence of randomness on thermodynamics is mainly rather typical. 
It leads to (i) a weak deformation of the entropy versus temperature curve 
with a decrease of the entropy at high temperatures (Figs. 4,5), (ii) a broad- 
ening and decreasing of the peak in the dependence of the specific heat versus 
temperature (Figs. 7,8), (iii) a smearing out of the cast in the mj versus Qq 
curve at T = for Qq = Jq and a nonsaturated transverse magnetization at 
any finite transverse field (Fig. 10), (iv) a decreasing and disappearing of 
the singularity (accompanying the saturation of fn^ at T — for Qq — Jq) in 
the curve Xzz versus f^o at T = (Fig. 12), and (v) a suppressing of static 
transverse linear susceptibility versus temperature curve (Fig. 13). 

However, as can be seen in Figs. 4-13, the influence of disorder, especially 
for small a, essentially depends on the sign of a. Particularly interesting is the 
case of strong asymmetry in the density of states p(E) when | a |Ri 1. ^From 
mathematical point of view the dependence of the computed quantities on 
temperature and averaged transverse field and the well-pronounced difference 
between the cases a ~ — 1 and a ~ 1 can be understood having in mind that 
these quantities according to (14)- (17) are integrals over E of products of 
p{E) (shown in Fig. 1) with functions with evident dependence on E at 
different f3. It is interesting to note that for some Hamiltonian parameters 
and temperatures even very large randomness (controlled by F) almost does 
not affect the thermodynamic quantities. This can be nicely seen in Figs. 
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4-13. 

It is worth to underline tiiat tiic asymmetry of p{E) leads to the ap- 
pearance of nonzero averaged transverse magnetization WT^ at zero averaged 
transverse field Qq. As it can be seen from (16) fn^ = at T = 0, fio = if 
dEp{E) = /q°° dEp{E). This is evidently true for a symmetric density of 
states p{E) (as in the case considered by H.Nishimori) but is not obvious in 
the considered case (13). The difference between the integrals j!^^dEp(E) 
and /o°° dEp{E) can be clearly demonstrated by numerical finite-chain calcu- 
lations as a difference between the numbers of negative and positive eigen- 
values Ap of the N X N matrix 1 1 Aij \ \ denoted by A/L and A/+, respectively. 
Examples for certain realization of the random model (l)-(3) are given in 
Tables I and II. The transverse magnetization for a certain realization at 
T = is given by = ^~2n'^ finds a good agreement of these 

direct numerical finite-chain calculations for —m^ with the results depicted 
in Fig. 10 (e.g. for T = 0.25 -ml ^ 0.051, 0.038, 0.011 if a = -1.01, -2, -5 
respectively, for T = 1 -w; ^ 0.095,0.030,0.003 if a = -1.01,-2,-5 re- 
spectively) . 

4 Conclusions 

In conclusion, we present exact calculations of the thermodynamics of the 
spin- 1 isotropic XY chain with random lorentzian intersite interaction and 

a transverse field that depends linearly on the surrounding intersite inter- 
actions (l)-(3). The derived exact expressions for the averaged density of 
states (13) and thermodynamic quantities (14)-(17) may serve as a testing 
ground for approximate methods of spin systems with off-diagonal disorder 
that usually involve an unclear error. Aside from this they are interesting in 
their own right, since experimentally accessible systems are always affected 
by randomness, and an understanding of disorder effects even within such 
simple model can help in comparing experimental observations and theoret- 
ical predictions. 

Unfortunately, the obtained results do not permit to calculate exactly the 
averaged spin correlation functions since such calculation requires the knowl- 
edge of averaged many-particle fermion Green functions. Spin correlations 
and their dynamics may be examined using exact finite-chain calculations 
developed in Refs. 26,27. 
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TABLE I. The numbers of negative and positive Ap, M- and A/V, for three 
reahzations of random intersite interactions in a finite model (l)-(3) [N — 
1000) with = 0, Jo = 1, r = 0.25. 





a = —5 


a = -2 


a= -1.01 


a = 1.01 


a = 2 


a = 5 




Ar_ 




Af- 




AA_ 


A/; 


Ar_ 


A/; 


Ar_ 




Ar_ 




0.997 


491 


509 


458 


542 


454 


546 


546 


454 


542 


458 


509 


491 


0.984 


490 


510 


464 


536 


442 


558 


558 


442 


536 


464 


510 


490 


1.008 


487 


513 


463 


537 


452 


548 


548 


452 


537 


463 


513 


487 



TABLE II. The numbers of negative and positive Ap, A/"- and A/+, for three 
realizations of random intersite interactions in a finite model (l)-(3) {N — 
1000) with Qo = 0, Jo = 1, r = 1. 





a = —5 


a = -2 


a= -1.01 


a = 1.01 


a = 2 


a = 5 




Ar_ 


A/; 


U- 




AA_ 




Ar_ 




Ar_ 




Ar_ 


A/; 


1.009 


495 


505 


470 


530 


402 


598 


598 


402 


530 


470 


505 


495 


0.986 


503 


497 


471 


529 


408 


592 


592 


408 


529 


471 


497 


503 


1.034 


494 


506 


469 


531 


406 


594 


594 


406 


531 


469 


506 


494 



List of figure captions 
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FIG. 1. The averaged density of states (13) p{E) vs. E — Qq. 

FIG. 2. The averaged density of states (13) p{E) (broken lines) and 
the density of states for certain reahzation of random intersite interactions 
obtained by exact finite-chain calculations (solid lines) vs. E — Qq. 

FIG. 3. The density of states p{E) averaged over 10 realizations for 
I a |< 1 obtained by exact finite-chain calculations vs. — Qq- 

FIG. 4. The entropy s (14) vs. temperature T = 0.25. 

FIG. 5. The entropy s (14) vs. temperature ^, F = 1. 

FIG. 6. The entropy s (14) vs. transverse field Qq, ^ = 0.1. 

FIG. 7. The specific heat c (15) vs. temperature ^, T = 0.25. 

FIG. 8. The specific heat c (15) vs. temperature ^, F = 1. 

FIG. 9. The specific heat c (15) vs. transverse field Qq, ^ — 0.1. 

FIG. 10. The transverse magnetization —rn^ (16) vs. transverse field 
at low temperature — 0.001). 

FIG. 11. The transverse magnetization — (16) vs. temperature ^ at 
Qo = 0.5. 

FIG. 12. The static transverse linear susceptibility —Xzz (17) vs. trans- 
verse field flo at low temperature (^ = 0.001). 

FIG. 13. The static transverse hnear susceptibihty — (17) vs. tem- 
perature 4 at Jlo = 0.5. 
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